Strongly Scale-dependent Non-Gaussianity 
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We discuss models of primordial density perturbations where the non-Gaussianity is strongly scale- 
dependent. In particular, the non-Gaussianity may have a sharp cut-off and be very suppressed on 
large cosmological scales, but sizeable on small scales. This may have an impact on probes of 
non-Gaussianity in the large-scale structure and in the cosmic microwave background radiation 
anisotropies. 
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Cosmological inflation [l| has become the dominant 
paradigm to understand the initial conditions for the 
Cosmic Microwave Background (CMB) anisotropies and 
Large Scale Structure (LSS) formation. This picture has 
recently received further spectacular confirmation by the 
Wilkinson Microwave Anisotropy Probe (WMAP) five 
year set of data Q. Present and future Q experiments 
may be sensitive to the non-linearities of the cosmological 
perturbations at the level of second- or higher-order per- 
turbation theory. The detection of these non-linearities 
through the non-Gaussianity (NG) [H has become one of 
the primary experimental targets. 

A possible source of NG could be primordial in ori- 
gin, being specific to a particular mechanism for the 
generation of the cosmological perturbations. This is 
what makes a positive detection of NG so relevant: it 
might help in discriminating among competing scenarios 
which otherwise might be indistinguishable. Indeed, var- 
ious models of inflation, firmly rooted in modern particle 
physics theory, predict a significant amount of primordial 
NG generated either during or immediately after infla- 
tion when the comoving curvature perturbation becomes 
constant on super-horizon scales |4j. While single-field 
[H and two(multi)-field 0] models of inflation generically 
predict a tiny level of NG, 'curvaton-type models' [7H9j, 
in which a significant contribution to the curvature per- 
turbation is generated after the end of slow-roll inflation 
by the perturbation in a field which has a negligible effect 
on inflation, may predict a high level of NG [lflj. Alterna- 
tives to the curvaton model are those models character- 
ized by the curvature perturbation being generated by an 
inhomogeneity in the decay rate ll|, [l2[ of the inflaton 
field. Other opportunities for generating the curvature 
perturbation occur at the end of inflation [3] and dur- 
ing preheating All these models generate a level of 
NG which is local, since the NG part of the primordial 
curvature perturbation is a local function of the Gaussian 
part generated on superhorizon scales. A phenomenolog- 
ical way of parametrizing the level of NG is to expand 
the fully non-linear primordial Bardeen gravitational po- 
tential, $, in powers of the linear gravitational potential, 



* = $L + /nl (#! - 
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The dimensionless quantity /nl sets the magnitude of the 
three-point correlation function 0] ■ In momentum space, 
the three point function (bispectrum), arising from the 
local NG is dominated by the so-called "squeezed" con- 
figuration, where one of the momenta is much smaller 
than the other two and it is parametrized by the non- 
linearity parameter J^l- Other models, such as DBI 
inflation 15 1 and ghost inflation [l6j], predict a differ- 
ent kind of primordial NG, called "equilateral", because 
the three-point function for this kind of NG is peaked 
on equilateral configurations, in which the lengths of the 
three wave- vectors forming a triangle in Fourier space are 
equal [l7| . The equilateral NG is parametrized by an am- 
plitude /nl" 1 ' 18 1- Present limits on NG are summarized 
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It is clear that detecting a significant amount of NG 
and its shape either from the CMB or from the LSS 
offers the possibility of opening a window into the dy- 
namics of the universe during the very first stages of its 
evolution and to understand what mechanism gave rise 
to the cosmological perturbations. Besides in the CMB 
anisotropies, where the Signal-to- Noise Ratio (SNR) for 
the NG is dominated by the small (angular) scales (3.[2l|. 
the primordial NG alters the clustering of dark matter 
halos inducing a scale-dependent bias on large scales [22j . 
Furthermore, non-Gaussianities are particularly relevant 
in the high-mass end of the power spectrum of perturba- 
tions, i.e. on the scale of galaxy clusters, since the effect 
of NG fluctuations becomes es pec ially visible on the tail 
of the probability distribution [23[. As a result, both the 
clustering properties and abundance of very massive ha- 
los are sensitive probes of primordial non-Gaussianities 
and could be detected or significantly constrained by the 
various planned large-scale galaxy surveys, both ground 
based (such as DES, PanSTARRS and LSST) and in 
space (such as EUCLID and ADEPT) [II]. 

The fact that the CMB anisotropies and the LSS probe 
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NG at different cosmological scales has recently spurred 
quite a lot of interest in the possibility that NG is run- 
ning with the cosmological scale 251- 27 1. For instance, 
single-field models with changing sound speed can have 
strongly scale-dependent NG [25j. Such models could 
evade the CMB constraints, but still have important ef- 
fects at scales responsible for the formation of cosmo- 
logical objects such as clusters and galaxies. A simple 
power law ansatz to describe the scale-dependence of 
NG is asumed for the bispectrum (<fr(ki)&(k 2 )$(k 3 )) = 
/nl (if/fcp)" /NL F(k l ,k 2 ,k 3 ), where K = (fcifc 2 fc 3 ) 1/3 , k p 
is some fixed pivot scale and n/ NL is the NG spectral 
index. Future observations such as CMBpol may be 
able to detect the scale-dependence of NG. The obser- 
vational error, around the fiducial value n/ NL = 0, is 
An /NL ~ 0.05(50// N L)/v/isk^ Hi) where / sky is the frac- 
tion of the observed sky. The reason why the power-law 
ansatz has been adopted in the literature is that the scale- 
dependence of NG turns out to be moderate in all mod- 
els analyzed, |fy NL | < 1. and usually proportional to the 
slow- roll parameters [27j . 

The goal of this short note is to show that the scale- 
dependence of the NG can naturally be much stronger 
than what have been discussed in the literature so far. 
The NG can be sizeable at small cosmological scales and 
totally negligible at large scales with the transition be- 
tween the two regimes taking place at some comoving 
scale ~ k^ 1 exiting the Hubble radius during the last 60 
(or so) e-folds of inflation. If one insists on making use 
of the NG spectral index, such a situation might be de- 
scribed by n/ NL ~ for comoving momenta k > fc* and 
n/ NL large and positive for momenta k < fc*. The rea- 
son why one is interested in such a possibility is readily 
explained. The square of the SNR corresponding to a 
non-vanishing NG bispectrum in the CMB anisotropics 
is proportional to 
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where, for simplicity, we have assumed the flat-sky ap- 
proximation, C(£) and B(£i,£ 2 ,£ 3 ) are the amplitudes 
of the two- and three-point temperature anisotropy cor- 
relators in multipole space, respectively. Consider, for 
example, the case in which the NG is dominated by 
the squeezed configuration (local case), l\ <C l 2 ,z- For 
£ ^> £-£, ~ 750 the combination of the radiation driving 
and Silk damping causes the exponential decrease of the 
angular power spectrum and bispectrum such that square 
of the SNR turns out to be proportional to 21 1 
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where £ m - m is a multipole of the order of £u and ^ ma x is 
the maximum multipole a given experiment can reach. 
Therefore, the SNR is dominated by the smallest (angu- 
lar) scale. The same is true for other shapes of NG; for 



instance for the equilateral NG, where the NG is a maxi- 
mized for £i ~ £ 2 ~ £ 3 , the square of the SNR scales like 
(/nl) 2 ^™- On the LSS side, the primordial local NG 
manifests itself by altering the clustering of dark matter 
peaks (halos) and inducing a scale-dependent Lagrangian 
bias between the power spectrum of halos and dark mat- 
ter, &l = Ph(k)/Pdm{k) on large scales [22| . The forecast 
from the shape of the large-scale power spectrum in the 
presence of NG in a given galaxy survey is sensitive to 
the largest scale ~ fc"^ within a given redshift slice and 
the square of the SNR scales like (f^l) 2 /k^ in ||. The 
LSS bias probe is sensitive to the large scale side of the 
primordial NG as the long wavelength part of the second- 
order perturbation modulates the position of the dark 
matter density peaks where galaxy reside. Therefore, a 
model of local primordial NG where the non-linearities 
are switched on at small cosmological scales and off at 
large scales will give no signal in the LSS bias, while it 
will possibly manifest itself in the CMB anisotropy. 

Let us now see how such a situation may be easily 
realized. Consider a scalar field o which gives a subdom- 
inant contribution to the energy density and which is ini- 
tially heavy with mass M ^> H, but which becomes light, 
m <C H at some (conformal) time r = r», corresponding 
to the time when the comoving wavelength 1 /fc* is exiting 
the Hubble radius. Here H indicates the approximately 
constant value of the Hubble rate during inflation. We 
will assume that the linear perturbations are dominated 
by the inflaton fluctuations, and since the er-field gives 
a sub-dominant contribution to the total energy density, 
the comoving curvature perturbation after the inflaton 
has decayed into radiation is ( « Crad+<W/ (4ptot), where 
we used £ = —H6p/p in spatially flat gauge, and that 
after inflation we have ptot = — 4i?ptot- Since the per- 
turbation in the radiation fluid seeded by the inflaton is 
approximately conserved, and can be computed the usual 
way during inflation, we only need to compute Sp a after 
the end of inflation in the two different cases: first for 
modes crossing outside the horizon when the mass of the 
cr-field is M H; secondly for modes crossing outside 
the horizon when the mass of the new field is m <C H. 

Let us first consider the modes which exited the hori- 
zon with mass M 3> H. The solution of the mode equa- 
tion during inflation is 



= -Ht 3 ' 2 
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(4) 



where v 2 M = 9/4-M 2 /iJ 2 (so that we have v M ~ iM/H 
for M > H) and m^'(x) are the Hankel functions. 
Picking up the Bunch-Davies vacuum for modes well in- 
side the Hubble radius gives C\ = (7r/2)exp (— irM/2H) 
and C 2 = 0. At late times, when (— kr) <C 1, the Hankel 
function behaves as 
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and therefore the solution for (— fcr*) <C 1 is 



2 M fc 3 I <zi? 



(6) 



After the mass becomes small the field almost freezes, 
and we have for super-Hubble scales at the end of infla- 
tion r e 3> T* 
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where i/ m ~ 3/2 — m 2 /(3H 2 ) and we will for simplic- 
ity assume that m 2 /(3H 2 ) > e = -H/H 2 . On the 
other hand, modes that exit the horizon after with 
(— fci"*) 3> 1 freeze in the moment they exit the Hub- 
ble radius in the usual way, and at the end of inflation 
on super-horizon scales they have the usual almost scale 
invariant spectrum given by 
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After the end of inflation when the inflaton has decayed 
into radiation, the Hubble rate will rapidly decrease, be- 
coming smaller than m, and the cr-field will quickly be- 
come non-relativistic. At this stage the amplitude |ok| 2 
will again begin to decrease like \/a . Thus, during radi- 
ation domination, for r 3> t , we can obtain the respec- 
tive solutions for |<7k| 2 averaged over an oscillation period 
by multiplying Eqs. and © by a factor (l/2)(a c /a) 3 . 

The vacuum expectation value of the cr-field will very 
quickly be driven to zero during the first stage of inflation 
when the mass of the field is large. Thus, the energy 
density fluctuations of the a field becomes 
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Note that the modes for which k <C a*H are sup- 
pressed by a factor H/M(k/a*H) 3 compared those for 
which k 3> a*H. However, even the modes that exit 
when the mass is small k 3> a*H give a small contribu- 
tion to the curvature perturbation at the end of infla- 
tion, but after the end of inflation they decrease slower 
than the background energy density of radiation and re- 
gain relative strength. Indeed, defining the power spec- 
trum (er 2 er 2 ) = (2ir) 3 6^(k + q)^V a 2(k) and setting 
rjcr = 2m 2 / (3H 2 ), we can write 

and following a calculation similar to that of Ref. [29j j . 
one can show that 
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which implies TV (ft) ~ (6/4r) a )(a c /a) 6 P 2 {k). 

The bispectrum of the curvature perturbation can be 
defined by 

(CfcXfeCfca) = (27r) 3 < 5( 3 )(k 1 +k 2 +k3)B c (fc 1 ,fc 2 ,fc 3 ) (12) 

such that a local type NG can be parametrized by a 
parameter by 



!?/nl C [P<(ki)Pdk2) + cyclic perm.] 



(13) 

where P ( = (27r 2 /fc 3 )P c . We will assume that the NG 
from the inflaton is negligible and that the fluctuations 
of the inflaton are uncorrelated with those of the cr-field. 
Under these assumptions, the bispectrum of the curva- 
ture perturbation is given by 
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In the squeezed limit fc 3 <C fci 



(15) 
it reduces to 
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The contribution to the curvature perturbation 
will be subdominant when the cr-field decays if 
(l/4)(6/256)(2e) 2 (a d /a ) 2 77 c7 7' 2 < V C , with the scale 
factor at the time of decay denoted by ad. Therefore, 
the cr-field maximally contributes to the spectrum say at 
the level of ~ 0(10)%, without disturbing the spectrum, 

if (a e /dd) ^ erfc^vV 2 ■ We then obtain that the NG is 
suppressed at large scales and is in the observationally 
interesting range 
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for comoving wavenumbers k > fc*. This simple set-up 
therefore predicts a NG which is basically vanishing at 
large scales and switched on at small cosmological scales. 
Let us now come back to the issue of how we can engineer 
a model where the mass of the cr-field changes abruptly 
at some time r* within the last 60 e-folds or so. Suppose 
that, beyond the inflaton and the a- field, there is a 
third field, \> such that the full potential is 
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V = V(<f>)+- (.gV + m 2 ) a 2 +- (-m 2 + h 2 ^ 2 ) X 2 + ^X 4 

(19) 
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If inflation is of the small-field type [l[ (it may be of 
the hybrid type or not), then in the first stage of in- 
flation the x ne ld acquires a vacuum expectation value 
(x) 2 = to 2 ,/ A, and the effective mass squared of the a- 
field becomes M 2 ~ (g 2 m 2 ,/X). It will be larger than 
H 2 if (g 2 /\) > (H/m x ) 2 . As soon as (<j)) 2 becomes 
larger than to 2 /h 2 the vacuum expectation value of the x 
field vanishes and the mass squared of the a becomes to 2 
which is assumed to be smaller than H 2 during inflation. 
The wavenumber fc* at which the NG are switched on 
corresponds to the scale exiting the Hubble radius when 
((f)) 2 ~ to 2 //i 2 . Notice that the same pattern can be ob- 
tained in large-field models of inflation (again, it might 
be of the hybrid type) , it suffices to have to 2 < and to 
change the sign of the 4> 2 x 2 interaction. 

Let us close with a couple of comments. First, the 
transition from the massive to the massless regime needs 
a fine-tuning to happen during the last 60 e-folds or so; 
nevertheless, if inflation is of the hybrid type, taking the 
parameters in the same ballpark of those added here, will 
naturally lead to a transition within the observationally 
interesting range of scales. Secondly, we have found it 
hard to build up models where the transition is from 
massless to massive regime, corresponding to a sizeable 
NG at large scales and suppressed at small scales: the 
NG perturbations generated on large scales tend to be 
overwhelmingly diluted in the subsequent evolution. Of 
course, one can turn this apparently negative comment 
upside down and conclude that strongly scale dependent 
NG will most preferably show up in the CMB. 
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